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Abstract 

We discuss the accuracy of the usual procedure for neutrino energy reconstruction which is based 
on the quasielastic kinematics. Our results are described in terms of a probability distribution for a 
real neutrino energy value. Several factors are responsible of the deviations from the reconstructed 
value. The main one is the multinucleon component of the neutrino interaction which in the case 
of Cherenkov detectors enters as a quasielastic cross section, increasing the mean neutrino energy 
which can differ appreciably from the reconstructed value. 

As an application we derive, for excess electron events attributed to the conversion of muon 
neutrinos, the true neutrino energy distribution based on the experimental one which is given in 
terms of the reconstructed value. The result is a reshaping effect. For MiniBooNE the low energy 
peak is suppressed and shifted at higher energies, which may influence the interpretation in terms 
of oscillation. For T2K at the Super Kamiokande far detector the reshaping translates into a 
narrowing of the energy distribution. 

PACS numbers: 13.15.+g, 25.30.Pt, 14.60.Pq 
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I. INTRODUCTION 



Neutrino oscillation experiments require the determination of the neutrino energy which 
enter the expression of the oscillation probability. This determination is commonly done 
through the charged current quasielastic events, commonly defined as those in which the 
emission product only include one charged lepton, the ejected nucleon being unobservable. 

For the "quasielastic" events where only the charged lepton is observed the only measur- 
able quantities are then its direction, i.e., its emission angle 9 with respect to the neutrino 
beam direction and its energy Ei (or kinetic energy T/ and momentum Pi). The neutrino 
energy Ey is unknown. The usual reconstruction procedure assumes that we are dealing with 
a genuine quasielastic event on a nucleon at rest. In this case the energy uj and momentum 
q transferred to the nuclear system are related by: u = {q^ — u"^) / (2M), where M is the nu- 
cleon mass. These quantities are related to the charged lepton observables by: uj = Ey — Ei 
and q"^ = Ey + Pf — 2EyPi cos 6. For illustration we take the example of an ejected muon. 
The quasielastic condition then gives the value Ey of the reconstructed energy 

E,-ml/{2M) 
' l-iE^-P,cose)/M- 

A binding correction can be introduced in this expression but it is irrelevant for our discus- 
sion. Several effects can influence this expression. Firstly the Fermi motion which broadens 
the quasielastic peak, the Pauli blocking which cuts part of the low momentum response, but 
also the fact that a number of events such as the multinucleon ones which are not genuine 
quasielastic ones can, in the case for instance of Cherenkov detectors, simulate those but 
have no reason to fulfill the quasielastic relation. The multinucleon component has been 



shown 



-5| to be responsible for a sizeable increase of the quasielastic cross section. 
In order to visualize the various effects on the reconstruction of the neutrino energy we 
first repeat an argument of Refs. 6] and j^. Eliminating E^ = uj + E^j^ in the expression of 
leads to the relation 

g2 _ ^2 ^ 4(^^^ ^ gi^2^ - + 2{E^ + uj){E^ - p^) cose. (2) 

In the u and q plane this equation represents a series of hyperbolas, which depend on the 
values of Efj_ and 6, with asymptotes parallel to the u = q line. Examples are shown in Fig. 
[H The intercept of the hyperbola with the region of response of the nucleus, whatever its 
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origin, quasielastic or not, fixes the possible u and q values for these values of and 9. 
In the case of a dilute Fermi gas where the region of response reduces to the quasielastic 
line, the intercept values uj = ^intercept and q = gintercept values are completely fixed. Hence 
the neutrino energy is also determined : Ey = + i^interccpt? which leads to the usual 
value of the reconstructed energy for this set of muon observables, E^ and 9. The Fermi 
motion introduces a broadening around the quasielastic line, as represented by the shaded 
area of Fig. [H We will show that in some situations this broadening is not innocent for 
the reconstruction problem. In addition the quasielastic peak itself can be distorted by 
Pauli correlations or by the collective nature of the response, as described in the random 
phase approximation (RPA) These effects destroy the symmetry of the distribution 
of the Ey values around the reconstructed energy E^. Moreover in a realistic description 
with correlations between nucleons, more than one nucleon can be ejected by the neutrino 
interaction. The Delta resonance can also be excited in this process with subsequent emission 
of several nucleons. This has the effect of enlarging the region of response to the whole uj 
and q plane. This means that for a given set of muon variables, E^j_ and 9, an infinity of 
UJ values are possible, hence of neutrino energies, instead of the unique quasielastic value 
implemented in the neutrino energy reconstruction formula. The question is to find the 
distribution of the possible E'^ values and to understand the influence of all these factors on 
the energy reconstruction. The results are described in terms of probability distributions to 
have a true neutrino energy between E^^ and Ey + dEy. Prior investigations by Benhar and 
Meloni [3] have dealt with the influence of the nuclear spectral function. Another work by 
Leitner and Mosel sl deals with the question of real pions produced and then absorbed on 
their way out of the nucleus, treated by transport equation. This process as well simulates 
a quasielastic interaction. In our work instead we omit the final state interaction, hence the 
absorption of real pions produced. 

Our work is based on the description of the neutrino cross section through the nuclear 



responses as is described in the work of Ref. [3| which successfully accounted for the Mini- 
BooNE data on the double differential cross section ^. The problem of finding the prob- 
ability distribution of the neutrino genuine energy in the "quasielastic" events cannot be 
dissociated from the neutrino energy distribution which is an input in this evaluation, as 
each energy is weighted with the neutrino flux, (^{Ey). In the applications of the present 
work we take three examples of the neutrino energy distribution, related to the correspond- 
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FIG. 1: (Color online) The neutrino hyperbolas defined by Eq. ([2]) for a muon kinetic energy 
T^=350 MeV, r^=750 MeV and several muon emission angles indicated in the figure. The shaded 
area delimited by the two corresponding continuous lines represents the region of the quasielastic 
response of a Fermi gas. The central dashed lines show the position of the quasielastic peak. 



ing neutrino oscillations experiments. One is for the MiniBooNE circumstances 10|, lllj . 
i.e. for this neutrino flux energy distribution that we take from Ref. jsj. For the second 
case we consider the T2K flux [12] at the near detector (ND), ^T2K{^pt)i which is somewhat 
more narrow. The energy profiles of the two beams are shown in Fig. [2J Finally in the 
application of our results to the energy distribution of the T2K electron events as observed 
at the Super Kamiokande far detector [l^, we take the electron neutrino flux ^T2K{^e) froin 
an evaluation in the oscillation description 



^T2Ki.^e) OC Sin^ 



Aml2L 



ND 
T2K 



(3) 



where |Am32| = 2.4 x 10 eV^ and L = 295 km. This electron neutrino flux is also shown 
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FIG. 2: (Color online) MiniBooNE and T2K normalized fluxes. The t'^ MiniBooNE flux is taken 
from Ref.jsl. The T2K flux at near detector is taken from [2]. The Ve T2K flux at Super 
Kamiokande far detector is obtained from Eq. ([3]). 



in Fig. |2l it is definitely narrower that the MiniBooNE one or the T2K near detector one. 



II. ANALYSIS AND RESULTS 



A. With specification of the lepton observables 

We first discuss what is the neutrino energy distribution f{E^,E^,9) for a given set of 
muon energy E'^ and angle 6. In this case one explores a single hyperbola in the {u, q) 
plane. The neutrino energy distribution follows that of the energy u transferred to the 
nuclear system along this hyperbola, as = + u, with a weight proportional to the 
neutrino flux at this energy ^{E^). Thus it is defined as 

HE,)dE,, (4) 

where C is a normalization constant. Several choices are possible for the normalization and 
we take here for C the inverse value of the total flux, C = {J dE^^{E^))~^, in such a way that 
the probability / has the dimensions of a cross section divided by an energy squared. The 



fiEjy, E^, 6)dE^ = C 



du dcos6 
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FIG. 3: (Color online) Probability distributions for T^=0.35 GeV evaluated at various values of the 
emission angle for the MiniBooNE flux. The vertical lines give the position of the corresponding 
reconstructed energy according to Eq. ([T]). 



neutrino cross section which enters Eq. (jlj) has two pieces, the one particle-one hole (Ip-lh) 
part, which is the genuine quasielastic part, and the multinucleon emission (np-nh) one which 
represents a sizeable fraction of the first one. The second part in particular introduces a high 
energy tail in the neutrino cross section, ignored in the usual reconstruction procedure which 
assumes that all the cross section for this energy and angle is concentrated at the intercept 
of the corresponding hyperbola with the quasielastic line. For all numerical evaluat^ns of 
this article we use the description of the nuclear responses of our previous work [3]. As 
in this work relativistic corrections to the nuclear responses are incorporated, the Ip-lh 
component of the response is treated in the random phase approximation (RPA) and we use 
the parametrization of the multinucleon component proposed in this work. In our evaluation 
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FIG. 4: (Color online) Probability distributions for T^=0.35 GeV evaluated at various values of 
the emission angle for the T2K near detector flux. The vertical lines give the position of the 
corresponding reconstructed energy Ei, according to Eq. ([1]). 

we do not introduce any experimental cut off which is maybe necessary to describe actual 
data. In this work we consider charged current neutrino cross sections on carbon which is 
the element constituting the MiniBooNE detector and the T2K near detector. In the case 
of Super Kamiokande the element involved is the oxygen but in Ref. ^j] we have shown that 
the cross sections per nucleon calculated for carbon and oxygen are almost identical both 
in the quasielastic and multinucleon channel. As a consequence we perform the calculations 
for carbon. 

The outcome of our study is illustrated in Figs. E] and H] which represent the probability 
distributions for one value of the muon energy, Efj_ = 0.45 GeV (T^ = 0.35 GeV) and 
evaluated at various values of the emission angle. Figure |5] represents the same quantities 
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FIG. 5: (Color online) Probability distributions for T^=0.75 GeV evaluated at various values of 
the emission angle for the MiniBooNE and T2K near detector fluxes. The vertical lines give the 
position of the corresponding reconstructed energy Ei, according to Eq. ([1]). 



for Tfj_ = 0.75 GeV. The corresponding hyperbolas for this particular muon energy and the 
selected values of the muon angle are displayed in Fig. [TJ We have performed the evaluations 
for the two neutrino fluxes, MiniBooNE and T2K at the near detector. The general features 
of the two distributions are rather similar in the two cases, with some differences that we 
will comment. They should be compared to the distributions for a dilute Fermi gas which is 
a delta function at the value of the reconstructed energy, E^, indicated by a line in Figs. |3l 
mandO The partial probabilities for Ip-lh and np-nh are shown separately, as well as the 
total one. The features of the probability distribution are radically different according to 
the value of the emission angle. The genuine quasielastic part, i.e., the Ip-lh part itself, can 
display sizeable deviations from the reconstructed energy. At small angles the corresponding 
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probability distribution is not symmetrical around E,y, and it is hardened. Part is due to 
the Pauli blocking effect which cuts the responses at small momenta. The RPA also has 
some hardening and quenching effect at small angles as it is displayed in Figs. |3l H] and [5] 
where the probability with and without RPA are shown. Moreover the np-nh part creates 
a high energy tail in the nuclear response, above the quasielastic peak, quite visible in the 
transverse part of the inclusive electron scattering data Larger energy transfer, u, 

meaning larger neutrino energies, the distribution acquires a high energy tail, absent for an 
uncorrelated Fermi gas, which produces an appreciable increase of the mean neutrino energy. 
With increasing muon angle the hardening of the distribution disappears at first and the 
Ip-lh distribution is centered near the reconstructed energy value. The np-nh component 
still creates the high energy tail. 

For larger angles an unexpected feature emerges, the energy distribution becomes softened 
with respect to the reconstructed energy. In the backward directions its distribution seems 
even unrelated to the reconstructed energy value. This feature can be understood from the 
corresponding hyperbola of Fig. [H With increasing muon angle the portion of the hyperbola 
inside the quasielastic region becomes very large. For instance for = 0.35 GeV already at 
a value of cos 6 = 0.2 the portion lies between the momentum transfer q ^ 0.6 GeV and g ~ 1 
GeV. In such a situation the form factors, vector or axial, have a strong suppression effect 
at the largest momenta (or energies), which explains the distortion of the neutrino energy 
distribution towards smaller energies which can be in addition favored by the flux factor. 
The problem becomes even more pronounced in the backward directions. For instance for 
cos^^ = —0.2 the distribution is concentrated at low energies and has little to do with the 
reconstructed energy, in particular for T2K (Fig. Hj). These large angles correspond to very 
large values of the momentum transfers, for which the np-nh component, which is difficult 
to evaluate reliably at these momenta, is not the main cause of the distortion effect. Notice 
that all the features mentioned above can be enhanced or diminished depending on the 
region of energy of the neutrino spectrum, which may or not favor the larger energies such 
as the np-nh tail. At small values, such as 0.45 GeV, the tail for a value cos 6 = 0.85 
occurs in an energy region of increasing flux for MiniBooNE and is favored while this is not 
so for T2K, which explains the difference in the importance of the tail. 
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FIG. 6: (Color online) Cosine of the muon angle solution of the Eq. ([5]) for various values of the 
reconstructed neutrino energy as a function of the muon energy. 

B. Distributions with no specification of the lepton observables 

In the previous study we have derived the probability distribution f{Ei,, E^, 0) for events 
with fixed muon variables: energy and emission angle. In the following we investigate a more 
global quantity, namely the neutrino energy distribution for a given neutrino reconstructed 
energy, F(E^, E^,), without any explicit reference to the muon variables. Many couples of 
values of these variables, muon energy and emission angle, can lead to the same reconstructed 
energy and one has to sum over these couples. A given value of the reconstructed energy. 
El, links the muon angle and its energy through 

KP^. cos e + M(e; - E^) - KEf, + ml/2 = 0. (5) 

We denote by cos 6{En, E,^) the cosine of the muon angle solution of this equation for a 
given set of values, E^^ and E„. Figure [6] displays its evolution with the muon energy for 
various values of the reconstructed energy. The probability energy distribution, F{E,^,Ei^), 
is obtained as an integral over the muon energy of the double differential neutrino cross 
section, taken at this particular value of the muon angle 
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FIG. 7: (Color online) Probability distribution before normalization, /, for several Ey values 
corresponding to the different vertical lines. The MiniBooNE and T2K near detector fluxes are 
used. 



(6) 

where the second identity defines the quantity J, i^min^max-j extreme values of i?^ 

obtained from Eq. ([5]) by making cos0 = ±1 and c is a normalization constant that we 
choose in order to have a total probability unity, / dEy F{Ey, E^) = 1, which gives to F the 
dimensions of an inverse energy. 

Some examples of these distributions are shown in Fig. [7] for some E^ values. Here 
we have plotted the quantity I{E,y,Ey) without the normalization factor c, instead of the 
probability density F, in order to single out the various contributions to the cross sections, 
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FIG. 8: (Color online) Probability distributions for several E,y values corresponding to the different 
vertical lines. Upper panel: using the MiniBooNE flux. Lower panel: using the T2K near detector 
flux. 



quasielastic with and without RPA, and the multinucleon contribution. The main features 
are similar to those found previously for fixed muon variables. This distribution has a high 
energy tail due to the np-nh cross section. For small or moderate E^, values the distribution is 
not symmetrical around E,,. The Ip-lh part itself is shifted towards larger neutrino energies 
by the Pauli blocking effects and the RPA influence. In addition the np-nh component is 
responsible for the high energy tail. For the smallest values of the reconstructed energy, 
^ 0.5 GeV, the np-nh tail is favored by the energy evolution of the neutrino flux which 
increases in this energy region, in particular for the T2K energy profile. With increasing 
values of the reconstructed energy the RPA hardening of the Ip-lh disappears and the form 
factor effect softens this part as explained previously. The np-nh compensates in part this 
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FIG. 9: (Color online) Probability distributions for several Ey values, corresponding to the different 
vertical lines, using the T2K flux predicted at Super Kamiokande for electron neutrinos. 



softening and the overall distribution looks approximately symmetric in particular with the 
MiniBooNE flux. For larger values E^, values the distribution becomes again asymmetrical 
but it is now the low energy side which is favored, largely due to the flux effect. These features 
are apparent in Fig. |H] where several probability distributions, F{Ey, E^), now normalized, 
are plotted in the same figure, showing the progressive evolutions with increasing values 
of the reconstructed energy. As illustrated in Fig. |9l this asymmetry effect is even more 
pronounced with the T2K far detector electron neutrino flux which is very narrow, as was 
shown in Fig. [2l The asymmetry of these distributions which favors higher energies at low 
El, values and smaller energies for large Ei, values has a narrowing effect for the smeared 
distribution as compared to the unsmeared one, as will be illustrated in the application 
to the T2K data. For the MiniBooNE case, with the assumption of a broad flux, the 
main effect is the transfer of low reconstructed energy events towards larger energies. This 
reshaping arises from the np-nh tail of the distribution. This is confirmed by the comparison 
of the normalized probability distributions with and without the multinucleon contribution 
as shown in Fig. [lOl For small values of the reconstructed energy the distribution in E^ with 
incorporates the np-nh contribution is sizeably deformed with respect to the one where this 
contribution is not included, it is suppressed at low E^, values and enhanced at large ones. 
This hardening effect is due to the important role of the np-nh tail in the normalization 
factor. It has interesting consequences for the analysis of the experimental data, as will be 
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FIG. 10: (Color online) Normalized probability distribution for some Ey values corresponding to 
the different vertical lines. The MiniBooNE flux and the T2K far detector Super Kamiokande z^e 
flux are used. 



discussed in the next section. 

The average neutrino energy, which depends on the reconstructed energy, ( -Ei/) average (-E'iy), 
is related to the normalized probability distribution F by 

( average = J dE„ E,y F{E^,E^). (7) 

The comparison with the reconstructed value gives an idea of the validity of the narrow 
quasielastic approximation. However in view of the shape of the distribution with the high 
energy tail, it would be meaningless to replace in the analysis of experimental data the 
reconstructed energy by the average one, only the energy distribution is significant. In the 
Table [T] some examples are given for the average neutrino energy of Eq. ([7]) for various 
values of E^ both for the MiniBooNE and the T2K near detector neutrino energy profiles. 

15 





MiniBooNE 




T2K ND 




(MeV) 
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300 


546 


335 


514 


350 


400 
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446 


500 


638 


527 


575 


528 
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800 


861 
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1000 
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981 


937 


914 


1200 


1190 


1164 


1116 


1104 



TABLE I: The average neutrino energy, Eq. ([7]), in MeV for various values of Ey considering the 
MiniBooNE and T2K near detector neutrino fluxes. The results are obtained with (QE+np-nh 
columns) and without (QE columns) the multinucleon contribution in the cross sections. 

We also give the same quantity but without np-nh contribution. Even in this case there is 
some slight modification at low values: the average neutrino energy is increased due to 
the RPA hardening of the responses. 

III. APPLICATIONS TO ACTUAL DATA ON ELECTRON NEUTRINO AP- 
PEARANCE 

Previous formulas given for muon emission apply for electron neutrinos, replacing by 
rrtf,. The results are in fact quite similar. For given lepton energy and angle the reconstructed 
energy is practically identical. For the same neutrino energy profiles the neutrino energy 
probability distributions are also quite similar in the two cases. In the following we apply 
our work to actual neutrino data of MiniBooNE jll| and T2K [13] on the observation of 
an electron excess in a beam of muon neutrinos. The experimental collaborations give the 
distribution of electron excess events in terms of the electron neutrino reconstructed energy. 
These neutrinos are attributed to the oscillation of muon neutrinos. As the oscillation 
phenomena depends on the neutrino energy it is important to have the distribution in terms 
of the real neutrino energy and not the reconstructed one which is the aim of this work. In 
a forthcoming one we will discuss the angular distribution problem and the distribution. 
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In our approach the transformation which leads to the true neutrino energy distribution is 
readily performed applying our previous probability distribution F. The number of events 
g{Ei,)dEi^ in a bin of reconstructed energy Ei, should be smeared by the function F{Ey, Ey) 
leading to the distribution G{Ey) in terms of the real neutrino energy 

G{E,) = j dE:g{E:)F{Ey,E:). (8) 

As the function F is normalized to unity this transformation conserves the area: 

j dKgiK) = J dE,G{E,). (9) 

We have performed this transformation both for the MiniBooNE and the T2K Super 
Kamiokande distributions. The evaluation of the smearing function F involves the electron 
neutrino energy spectrum. In the T2K case we have given the evaluation of this quantity 
from the muon neutrinos spectrum and the oscillation formula of Eq. ([3]). For MiniBooNE 
we make for the moment the assumption that the energy distribution is the same as the initial 
muon one. In the future we will improve this approximation in an oscillation model. As 
a last remark notice that our procedure is reversible: starting from a neutrino distribution 
in terms of the true energy, for instance a theoretical one as will be given for T2K, we 
could apply a similar procedure to derive the corresponding distribution in terms of the 
reconstructed energy, directly comparable to the experimental data. 



A. MiniBooNE 



We first start with the MiniBooNE case where the so called "low energy MiniBooNE 
anomaly" has been displayed Q, with a number of electron event excess concentrated 
in the low energy region of the reconstructed energies. The MiniBooNE collaboration [l^ 
has proposed to discard the low energy points responsible for the anomaly, which do not fit 
easily the oscillation interpretation. But there is no obvious reason to do so and here we 
keep all the points which possibly represent an electron excess. The experimental results are 
given in terms of an excess number in large intervals of E^ (200-300MeV, 300-375 MeV, etc). 
In view of this and for simplification, instead of the continuous integral on Ey of Eq. (|8]), 
we have taken for F a discrete sum of three energy points in each bin, the extreme and the 
central ones (which explains the spikes of our curves). The result of our study is contained in 
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FIG. 11: (Color online) Three different treatments of experimental MiniBooNE data 
neutrino excess events. The histogram is obtained identifying the reconstructed neutrino energy 
Ely with the true neutrino energy Eu. The full line is our smeared distribution including the 
multinucleon contribution. The dashed line is our smeared distribution without the multinucleon 
contribution. 



Fig. [TTJ The full curve represents our smeared distribution obtained from the experimental 
data excess [ll|. The experimental histogram is also shown in Fig. [TT] under the assumption 
Ejj = Ejj. For information we also give the smeared distribution in the absence of the np- 
nh contribution (dashed curve) which is more similar to the original data. However for us 
the more realistic curve is the one with the inclusion of the np-nh component. We stress 
that all the distributions shown in Fig. [11] represent the same experimental data with the 
different treatment of the abscissa, the neutrino energy. The areas under these curves are 
the same. The salient feature is that our smeared distribution (with the np-nh component) 
is broadened and hardened, with the suppression of the low energy excess, as compared to 
the initial histogram. The effect arises from the high energy tail of our smearing function. 
The difference with the original distribution is so large that it could affect the interpretation 
in terms of oscillation. Our smeared distribution with the np-nh component is closer to 



the LSND 



15j solutions shown in Ref. ll| than the original histogram. At present this 



important conclusion is only tentative as the flux evaluation should and will be improved. 
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1. Discussion and tests 



In view of the possible implications of this work for the neutrino oscillation interpretation, 
several questions have to be raised. We have pointed out the uncertainties of the present 
evaluation of the smearing function linked to the choice of the electron neutrinos energy 
distribution. This will be improved. Another point is that in the MiniBooNE experiment 
there is a large experimental error on the number of interesting events linked to the existence 
of a large background which has to be subtracted. This subtraction is the reason for the 
negative value of the distribution in some bins, that we have kept in our evaluation. This 
background is also given as a function of the reconstructed energy Ey. What part of the 
background should be smeared in the same way and what would be the influence on the 
distribution of interest? This should be for instance the case for the background due to 
electron neutrinos originally present in the beam. This question should be cleared. 

On the theoretical side the main question is : how reliable is our smearing ? The hardening 
effect arises from the high energy tail, due to the multinucleon contribution, of the smearing 
function for low values. It may change somewhat with different energy distributions of 
the electron neutrinos but we believe that the suppression of the low energy peak with a shift 
towards higher neutrino energies is a robust prediction that will survive the better description 
of this distribution. Independently of the flux problem how reliable is our evaluation of this 
part? We remind that the basic input for the evaluation of the smearing functions is the 
double differential cross section, -^^^^f^frg- We have tested this quantity jsj] in connection with 
the experimental MiniBooNE data on the double differential cross section with respect 
to the muon observables, where it also enters. Our theoretical evaluation of '^^dcoso 
been quite successful in the reproduction of these data. The multinucleon component is an 
important part of this cross section and is needed to reproduce the data in all the kinematical 

n n 

domain [5[,p|]. This double differential cross section involves a different weighting than the 
distribution F, by the muon neutrino MiniBooNE energy distribution 



d'^a 1 



dEf, dcos9 f <^{Ey) dE^ 



dE,. 



d^a 



dco dcosO 



HE,). (10) 



Hence for a given value of E^ several energy transfer are involved. Although the link with the 
present evaluation is therefore not a direct one we can argue that the success of our previous 
description {sl can give confidence in the present evaluation, as explained in the following. 
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0.2 GeV < T < 0.3 GeV 0.3 GeV < T < 0.4 GeV 




cos 9 cos G 

FIG. 12: (Color online) MiniBooNE flux-averaged charged current "quasielastic" f^-^^C double 
differential cross section per neutron for some values of muon kinetic energy as a function of the 
scattering angle. The theoretical pure quasielastic (dashed curve) and the (solid curve) generalized 
quasielastic (QE+np-nh) cross sections are those obtained in The MiniBooNE experimental 
points are taken from Q]. 



The suppression of the low energy peak of the MiniBooNE distribution of Fig. [TT]comes from 
the behavior of the smearing function for low values of the reconstructed energy {E,, ~ 0.2 
to 0.4 GeV) with its high energy tail (see the corresponding distributions of Figs. IHl and fTOll . 
We thus consider the example of the smearing function for the value of the reconstructed 
energy, E^, = 0.3 GeV, and a large value of the real energy, E^, =0.7 GeV, which is well in 
the tail region and hence dominated by the multinucleon component that we want to test. 
We investigate what is the relevant range of the energy transfer for these two values of the 
neutrino energies. The bounds for the electron energy are visible in Fig. [6] in which the 
abscissa is now the electron energy. They correspond to the extreme values of cos^ = ±1. 
Neglecting the electron mass the bounds for its energy E^ are E^/{1 + {2E^/M)) < E^ < E^, 
or 0.2 GeV< Ee < 0.3 GeV, which leads for the bounds of the energy transfer, 0.4 GeV< 
u = {El, — Ee) < 0.5 GeV. The question is then if this region is covered by the data on the 
double differential cross section -tb^ — ^ measured by MiniBooNE. In this cross section as 
given in Eq. ( fTOl) . the muon neutrino flux, peaked at E,y ~ 0.7 GeV, favors the region of 
muon energy in the vicinity of E^j, = 0.7 GeV — u. For 0.4 GeV< u <0.5 GeV the region 
favored by the peak position is between i?^ = 0.2 and 0.3 GeV. As the muon neutrino energy 
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peak is broad, somewhat larger values are relevant as well to test the region of energy 
transfers between 0.4 and 0.5 GeV. We select therefore for the test of our description of the 
np-nh component in the region of interest the two smallest values for experimentally 
available, E^ =0.35 GeV and E^ =0.45 GeV. The comparison with the MiniBooNE data is 
shown in Fig. [T^l reproduced from our previous work [3]. The np-nh part is clearly needed 
to account for the data. 

We can therefore conclude that there are experimental indications to confirm the impor- 
tance of the high energy np-nh tail and hence the reshaping towards higher energies of the 
MiniBooNE distribution of the low energy electron excess when expressed in terms of the 
real neutrino energy rather than of the reconstructed value. 

B. T2K 

The electron neutrino energy distribution at the far Super Kamiokande detector as cal- 
culated from the oscillation expression has been displayed in Fig. [21 The spectrum is much 
narrower than the presumed one for the MiniBooNE electron excess events. Hence the re- 
shaping features are somewhat different. On the experimental side the detection of excess 
electron is still in an early stage and some candidates (6) have been detected. The six elec- 
tron events are spread on large bins of reconstructed energy. Their energy distribution is 
shown in the histogram of Fig. [T3] (assuming Ei, = E^). We also show in this figure our 
smeared distribution with the same methods as for MiniBooNE. The main feature here is 
the narrowing of our energy distribution as compared to the experimental histogram which 
assumes the identity of the reconstructed and the true energy. We have explained the origin 
of this narrowing from the different behavior of the smearing function at small and large 
El, values. In the same figure we also show our theoretical prediction for this distribution. 
Since it represents the interaction number for the electron neutrinos which underwent a 
"quasielastic" event, it should be given by the quantity crQE+np-nh(-E'i/J$ff)^ (-Ei^J, product 
of the total "quasielastic" cross section by the electron neutrino energy distribution. The 
area below the theoretical curve has been normalized to the same total number of events. 
We take for the cross section, crQE+np-nh(-Ei/g), our calculated value, similar to the one for 
muon neutrinos [1,] which accounted quite well for the MiniBooNE data. Notice that, due to 
the normalization to the number of electron events, the introduction of the np-nh component 
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FIG. 13: (Color online) The experimental T2K data on electron neutrino excess events. The 
histogram is obtained identifying the reconstructed neutrino energy E,y with the true neutrino 
energy E^- The full line is our smeared distribution including the multinucleon part. The dot- 
dashed line is our theoretical prediction, i.e., the product aQE+np-nh{Eue)^j^2Ki'^i'e) with an area 
normalized to the total number of events. 

in the "quasielastic" cross section has practically no effect, since it amounts essentially to 
an increase of the cross section by a constant multiplying factor. It is in the narrowing of 
the smeared distribution that the np-nh shows up. In spite of the small number of events 
it appears that the narrowing goes well in the direction of the theoretical prediction. This 
offers a consistency test. If the electrons were background ones there would be no reason 
for an agreement. A similar comparison will be performed for MiniBooNE when a better 
smearing function will be evaluated. 

IV. CONCLUSION 

We have examined in this work the validity of the approximation contained in the identi- 
fication of the reconstructed neutrino energy with the real energy. The errors introduced by 
this assumption have various origins. The straightforward ones are the Pauli blocking effect, 
the Fermi motion which broadens the region of response and the collective aspects of the 
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response. However the main correction arises from the existence of the multinucleon com- 
ponent in the neutrino cross section on nuclei which can in Cherenkov detectors simulate a 
quasielastic process. We have introduced the probability distribution to have a real neutrino 
energy E^, given a reconstructed energy value E,^. We have found that at low Ey values this 
distribution acquires an important high energy tail due to the multinucleon component. For 
large reconstructed energies it is instead the low energy tail which is favored. The neutrino 
energy distribution is an important ingredient in these evaluations since it appears as a 
weight factor in the cross sections. Its influence can be concisely described as follows: the 
transition from the reconstructed to the real energy tends to concentrate the events in the 
region of high flux provided there is some strength in this region. 

For the MiniBooNE data on the energy distribution of the excess electron events, with 
some simplified assumptions on the electron neutrinos energy distribution, we find that it is 
modified with a suppression of the low energy peak, shifted at higher energies. This feature, 
when confirmed, will have consequences for the compatibility with other short-baseline data, 
in particular with the LSND ones. In this respect one question will have to be answered: 
if other short-baseline experiments should not be affected in a similar way. There is no 
general answer and the type of neutrino reaction as well as the detection method should be 
examined case by case. 

A reshaping occurs as well for the T2K electron events in the Super Kamiokande 
detector. In this case the narrowness of the electron neutrino energy distribution shows 
up as a narrowing of the electron events distribution, which goes in the direction of our 
theoretical prediction. In a forthcoming work we will also investigate the corresponding 
problems for the muonic events both in the close and in the far detector. In the latter 
case the oscillations lead to an energy distribution of the muon neutrino spectrum which is 
quite different from the one at the close detector {l2]. The reshaping in terms of the true 
neutrino energy may then reveal some amusing and unexpected features. But we also keep 
in mind the problems associated with the MiniBooNE data and what we can say at the 
nuclear light on the challenging problems of sterile neutrinos and short-baseline oscillations. 
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